Optimal operation of transition-edge sensors on ballistic membranes 
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We calculate the operating parameters of a transition edge sensor that is mounted on a thin dielectric 
membrane with the assumption that the phononic heat transport in the membrane is ballistic. Our 
treatment uses the correct phonon modes from elasticity theory (Lamb- modes), and spans the tran- 
sition from 3D to 2D behavior. The phonon cooling power and conductance have a global minimum 
as function of membrane thickness, which leads to an optimal value for the membrane thickness 
with respect to noise equivalent power at a fixed operating temperature. The energy resolution 
of a calorimeter will not be affected strongly, but, somewhat counterintuitively, the effective time 
constant can be reduced by decreasing the membrane thickness in the 2D limit. 
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Superconducting transition edge sensors (TES) are 
currently under heavy development to be used as ultra- 
sensitive calorimeters and bolometers over a wide range 
of frequencies, from sub-mm radiation to gamma raysii 
Some recent examples of ambitious projects using TES 
sensors are: the new sub-mm camera for the James Clerk 
Maxwell telescope (SCUBA-II)^, the detection of sin- 
gle near-infrared photons for quantum cryptography^, 
the X-ray imaging spectrometers for future ESA, NASA 
and JAXA mission*!, and the 7-ray detection of nu- 
clear materials^. In most of these detector designs, the 
superconducting TES film has been thermally isolated 
from the surroundings by mounting it on a thin dielec- 
tric membrane, usually made of amorphous silicon nitride 
(SiNx) due to ease of fabrication. This membrane lim- 
its the thermal conductance to the bath, and is therefore 
critical for the operation of the devices. 

The phonon transport in the membrane can be either 
diffusive or ballistic, depending on how easily phonons 
are scattered in the sample. At the low temperatures 
where TESes are operated (^ 0.1 K), the bulk scatter- 
ing mechanisms (mass impurities, phonon-phonon scat- 
tering, phonon-two-level system scattering) become very 
weak)^ leading to surface limited thermal conduction. If, 
in addition, the surfaces are smooth on the length scale 
of the dominant thermal phonons. the surface scatter- 
ing is mostly specular, and phonon transport becomes 
ballisticii Ballistic phonon transport has been observed 
for crystalline bulk samples a long time ago^, but was also 
recently shown to be valid for some SiNx membranes2ii£. 
For thermal sources such as TES film radiators in the 
ballistic limit, the analogy with photon thermal black- 
body radiation is apparent if the substrate is three- 
dimensional. Then the emitted power has the typical 
Stcfan-Bolzmann form P = AaT'^, where for phonons 
a = TT^ kg / {15h^)'Eei / cf , summing over the different 
phonon modes with speeds of sound Ci and radiator emis- 
sivities e^J^ On the other hand, it is less clear what hap- 
pens to ballistic phonon transport when the substrate 
becomes two-dimensional. This is a realistic concern, as 
the dominant emitted phonon wavelength is of the order 
of 1/Ltm for SiN at 100 mKfi^ comparable with a typical 
membrane thickness, so that many practical devices are 



at least approaching the 2D limit. 

In this paper, we describe a theory for operating TES 
detectors on thin membranes, spanning the transition 
from a fully 3D to a fully 2D substrate (membrane), us- 
ing elasticity theory. The eigenmodes of a thin mem- 
brane are no longer the usual plane wave phonons, but 
are so called Lamb modes, with non-trivial displacement 
fields and dispersion relationsiii For this reason, ther- 
mal conductions^ and the phononic heat capacityi^ are 
strongly affected, leading to the interesting effect that 
radiated power and thermal conductance have an abso- 
lute minimum as a function of the membrane thickness. 
Because of this effect, the noise equivalent power (NEP) 
also has a minimum, leading to the notion of an opti- 
mal membrane thickness for TES bolometers. However, 
the effective time constant also has a maximum at the 
same point, thus, one needs to make sure that the nega- 
tive effect of slowing the detector down is not critical for 
the application. We also calculated the influence of thin 
membranes on the phonon noise limited energy resolution 
of TES calorimeters, and found that the effect is weak, 
especially around the optimal detector temperature. 

In isotropic 3D bulk systems there are three indepen- 
dent phonon modes, two transversally and one longitu- 
dinally polarized, with sound velocities Ct and ci , respec- 
tively. In the presence of boundaries, the bulk phonon 
modes couple to each other and form a new set of eigen- 
modes, which in the case of a free standing membrane 
are horizontal shear modes (h) and symmetric (s) and 
antisymmetric (a) Lamb modes^^. The frequency uj for 

the h modes is simply uj = Ct^k"^^ -\- {mn/d)'^, where k\\ 

is wave vector component parallel to membrane surfaces, 
d is the membrane thickness and the integer m is the 
branch number. However, the dispersion relations of the 
s and a Lamb modes cannot be given in a closed analyt- 
ical form, but have to be calculated numerically^^ The 
lowest three branches, dominant for thin membranes at 
low temperatures, have low frequency analytical expres- 
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ity of the s mode, and m* = Ti 

an effective mass for the a-mode "particle". This lowest 
a-mode with its quadratic dispersion is mostly responsi- 
ble for the non-trivial behavior for the detector perfor- 
mance described below. 

To simplify the discussion, wc assume that the ther- 
mal conductance is only limited by the membrane itself, 
in other words we do not consider the effects of electron- 
phonon non-equilibrium or thermal gradients within the 
TES filmic, or of boundary resistance between the TES 
film and the membrane, and the membrane and the sup- 
porting 3D substrate. The importance of these added 
effects depends critically on the materials and detector 
geometry, and can in principle be minimized. With these 
assumptions the total heat flow out of the detector is 
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FIG. 1: (a) The differential thermal conductance g of a bal- 
listic radiator as function of the membrane thickness d. The 
crossover thickness dc = hct/2kBT. g is normalized with 
po = {kslct /8Tv^){kBT /hctf . The asymptotic 2D and 3D 
limits, computed from Eqs. [3] and |4] are indicated by the 
dashed lines, (b) The net emitted power P„et of a detector 
at temperature Td connected to the substrate at temperature 
Ts as function of membrane thickness d, for temperature ra- 
tios Td/TS = 1.01, 1.1, 1.2, 1.4 and 2.0. Curves with higher 
Td/Ts lie above curves with lower Td/Ts- ds = hct/2kBTs 
and Po = l(kBTsf/{27vh^ct). 



where I is the circumference of the detector, n{uj, T) is the 
Bose-Einstein distribution and a and m are the mode and 
branch indices This expression can be used to compute 
the transition from 3D to 2D if enough branches arc used. 
If the membrane is thin and temperature low, i.e. Td <C 
hct/2kB, only the lowest branches (Eqs.[T|) are occupied, 
and we are fully in the 2D limit, in which case 
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Note that the effective mass of the lowest a mode depends 
on the membrane thickness and hence in the 2D limit 
P oc 1/Vd. In the 3D limit Td ^ hct/2kB, the dominant 
phonon wavelength is much smaller than d, leading to 
decoupling of the longitudinal and transversal modes and 
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As expected, P3D oc d. This means that decreasing the 
membrane thickness at a fixed temperature, the radiated 
power will first decrease with d, then reach a global min- 
imum and will increase again, if we decrease d further. 
The minimum is approximately at the 2D-3D crossover 
thickness dc = Tict/ {2kBT). The same behavior occurs 
for the differential thermal conductance g = dP/dT. In 
the 2D-3D crossover range, we have computed P and g 



numerically, using the lowest 100 branches. Fig. [IJa), we 
plot g as a function of d, showing an increase of 5 by ~ 
factor of 5 when decreasing d from dc to IQ-^dc, cor- 
responding to changing d from 240 nm to 2.4 nm at 100 
mK for SiN^. 

The net phonon cooling power of the detector is the dif- 
ference between the power flow out of the detector and 
the power flow into the detector that is emitted from the 
substrate. Hence, if we denote the detector temperature 
by Tr) and the substrate temperature by T5, the total 
cooling power is Pnct PiTo) ~ PiTs)- In Fig. W^h) we 
plot Pnct as function of d for different ratios To/Ts- As 
we have two different temperatures, we chose to normal- 
ize the membrane thickness with respect to Ts and define 
the quantity ds = hct/2kBTs- The thickness dependence 
is the same as for g, and as one would expect, the radi- 
ated power increases with increasing Tp. The place of the 
minimum is slightly shifted from curve to curve, as the 
2D-3D crossover is dominated by the detector phonons 
at high Td/Ts. 

A useful way to operate TES detectors is by voltage 
biasing. This leads to negative electrothermal feedback 
(ETF), which speeds up the detectors and improves their 
sensitivity and energy resolution^. For TES on thin bal- 
listic membranes, the basic operational theory applies, 
as long as one takes into account the correct formulas for 
Pnet and g as discussed above. Especially interesting is 
the question: Will the non-trivial thickness dependence 
lead to observable effects for the important detector pa- 
rameters such as the noise equivalent power NEP, effec- 
tive time constant TeS and energy resolution AE7 In 
the following we discuss this issue, noting that we only 
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consider here the contributions from the thermodynamic 
noise generated by the phonon thermal transport in the 
membrane, the so called phonon noise. Other noise mech- 
anisms can be easily added to the discussion using known 
formulas^. 
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FIG. 2: (a)The ballistic phonon NEP of a TES as func- 
tion of d for temperature ratios Td/Ts = 1.0, 1.1, 1.2, 1.4 
and 2.0, ordered from bottom to top. The normalization 
constant is NEPo = ^/l/c't{kBTsf/{2nh). (b) The effec- 
tive time constant Toff as function of d for the same ratios 
Td/Ts = 1 (solid line) and 1.1, 1.2 and 1.5 for a = 10 
(dashed lines) and a = 100 (dotted lines). Curves with higher 
Td/Ts below curves with lower Td/Ts for the same value of 
a. reff,o = 8TT^h^ct'yV,i/ilk%Ts). 

The phonon noise power spectral density for bal- 
listic transport can be calculated from microscopic 
considerations^SiSi leading to an equation for the phonon 
NEP: 

NEP = ^2kB[Tlg{TD) + TlgiTs)]. (5) 

In Fig.[2^a), we plot the NEP as a function of d. It has a 
clear minimum, just like g, as expected from the simplic- 
ity of Eq. . In the 3D limit the NEP becomes propor- 
tional to Vd, whereas in the 2D limit it approaches d~^^'^. 
The absolute minimum in the NEP means that it is possi- 
ble to define the optimal membrane thickness dopt, which 
gives the highest sensitivity for a TES bolometer at fixed 
Tjj and Tg. For an equilibrium bolometer (To = Ts) 
dopt ~ ds, whereas for Td > Ts, dopt ~ hct / {2kBTD) ■ 
For SiNx membranes and To = Ts = 0.1 A', dopt = 240 
nm. 

In addition to sensitivity, an important characteristic 
of a bolometer is its response time Tgs- For a volt- 
age biased TES of volume V with its electronic heat 
capacity Cei = jVTd being dominant, in the perfect 
voltage bias and low inductance limil^i it is given by 
Teff = Cei/{g{TD) + aPnet/To), where a = T/RdR/dT 
is the steepness parameter of the transition under bias. 
Fig. [2ljb) presents TeS as function of d for several val- 
ues of Td/Ts and two examples of a = 10,100. Now 
Teff has a maximum around hct/{2kBTu) and is propor- 
tional to -y/d in the 2D limit and to l/d in the 3D limit. 



Thus, if one wants to operate near dopt, limitations of 
the time constant need to be considered. Fortunately, 
biasing strongly into the ETF (high a), helps to reduce 
Toff significantly. 

In Fig. [3l^a) . we show the NEP as function of Td for 
different values of d. In general the NEP increases with 
To- For a fixed value of To/Ts we can always find two 
membrane thicknesses that yield the same NEP, one with 
d < ds and one with d > ds- However, the To depen- 
dence of these two values are quite different, as can be 
seen by comparing the slopes of the curves at crossing 
points. 
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FIG. 3: (a) The phonon NEP vs. Td/Ts for membrane thick- 
nesses d/ds = 0.004, 0.04, 0.4, 1, 4 and 40. The curves for 
d/ds < 1 are plotted with solid lines with higher curves corre- 
sponding to lower d, while dotted lines represent curves with 
d/ds > 1 with lower curves corresponding to lower d. (b) 
The high C FWHM energy resolution AEfwhm of a TES 
as function of the ratio Td/Ts for the membrane thicknesses 
d/ds = 0.04, 0.4, 1, 4 and 40, ordered from the bottom to the 
top. AEo = 2-8y/'yVkBTl/ai- 

Finally, the phonon noise limited full width half 
maximum (FWHM) energy resolution of an optimally 
filtered TES calorimeter— for high loop gain C ~ 
aPnet/i9iTo)To) » 1 is 

I 2 / NEP'^ \ 

A£;fwhm = 2.355JkBTlCei^ ^ „ , 

V ^ ^ \4kBToPnetJ 

(6) 

where aj = T/RdR/dT- AiJpwHM is plotted as a 
function of To in Fig. El^b) for different membrane thick- 
nesses. It has a minimum just below To/Ts ~ 1.2, al- 
most independent of d and in agreement with the 3D 
ballistic result^^, so that an optimal To for a fixed Ts 
can be defined. For higher ratios To/Ts, the thinner 
membranes decrease AiJpwHM slightly, by about 6 % at 
To/Ts = 2 from d = 40ds to d = O.OAds- 
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